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Abstract 

The CDF 3.7 sigma evidence of pure annihilation B s — > ir + ir~ decays, if confirmed, would imply 
a large annihilation scenario in the QCD factorization approach. This is somewhat unexpected as 
the large annihilation scenario was disfavored in previous studies. In this paper we reinvestigate 
the role of annihilation topology in QCD factorization. We find that it is not easy to reach the 
CDF central value of B s — > 7r + 7r~ decays when other decay channels are considered. Our analysis 
also reveals that, for well- measured charmless B decays into two final pseudoscalar mesons, the 
QCD factorization predictions with large annihilation parameters show good agreement with the 
experimental data except B s — > K + K~ and — > K°K° decays. Though other possibilities 
can not be excluded, this may indicate that the SU(3) flavor symmetry breaking should be taken 
into account for the annihilation topology. In addition, there are different annihilation topologies, 
so that somewhat different annihilation parameters may be chosen for different final states and 
different annihilation topologies. If so, the predictive power of the QCD factorization method may 
be rather limited for many decay channels. 
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Very recently the CDF collaboration has reported a 3.7 a evidence for pure annihilation 
B s —7- 7r + 7r~ decays, together with a measurement of (also pure annihilation) Bd — > K + K~ 
decays []| 

B(B S ->• 7T+7T-) = (0.57 ± 0.15 ± 0.10) x 10~ 6 , 
B{B d -> K + K~) = (0.23 ± 0.10 ± 0.10) x 10~ 6 , (1) 

where the first errors are statistical and the second systematic. As the large hadron collider 
(LHC) has been running very well this year, the LHCb collaboration should be able to check 
and improve these results very soon. 

Theoretically it was noticed first in {2- 4] that weak annihilation amplitudes may not 
be negligibly small in charmless B decays and was predicted in perturbative QCD method 
111 with the same central value as the experiment. In QCD factorization method 

(QCDF) 8l4ll|. although weak annihilation contributions are formally power suppressed in 
A.Q C D/ m bi they are supposed to be important, together with the chirally-enhanced power 
corrections, to account for the large branching ratios of penguin-dominated B decays. In 
addition, the annihilation topologies may provide large strong phases which are crucial to 
accommodate the significant direct CP violation of Bd — > tt~K + . Unfortunately these 
power correction terms are not calculable in QCDF as they contain endpoint singularities 
which violate the factorization theorem. Phenomenologically these chirally-enhanced power 
corrections and weak annihilation amplitudes are estimated in a model dependent way 10] , 



and comprehensive studies of charmless B decays (see, e.g., llMl6l|) show good agreement 



with the experiment in general. However in soft collinear effective theory (SCET) 17H19|. 



it was argued that the chirally enhanced power corrections and weak annihilation diagrams 



are actually factorizable and real 



20 



2l| , while the so-called charming penguins 22j-[24j are 



supposed to be non-perturbative and important to account for the large branching ratios 
and CP violations of charmless B decays. 

Noticed that theoretically it is still controversial on whether the charming penguins would 



invalidate the standard picture of QCD factorization [25|, |26(. It is also hard in practice to 
tell whether charming penguins are really important, as generally weak annihilation and 
charming penguins have the same topology Q. But with the first experimental evidence of 
pure annihilation B s — > 7r + 7r~ decays, it is clear now that the annihilation contributions 



1 The ratio B(B S — > p + K )/B(B s — > n + K ) may provide some insight into this issue, as discussed in 



27|. 
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must be significant in charmless B decays, irrespective of whether the charming penguins 
are large or not. 

It is therefore timely to reinvestigate the annihilation contributions in QCDF, especially 
considering that the annihilation amplitudes seemed to be underestimated in previous QCDF 
studies. For instance, B(B S — > 7t + tt~) was estimated to be 0.155xl0 -6 in the favored scenario 
54 of ll(|, which is about three times smaller than the CDF observation. One might wonder 
that, with the annihilation magnitude larger than expected, whether it is still possible for 
the QCDF predictions of charmless B decays to be consistent with the experiments. This 
motivates our work below. 

In the framework of QCDF, the decay amplitudes of B s — > and B d — > K + K~ can 

be expressed as TlJ 



A(B d -> K + K~ 



Ktt [ Kb V us 



B KK ( V ub V ud 



hi + 26 4 + ^b 4 ,EW 

b x + 26 4 + ~b 4 ,Ew 



+ VSy„ 
+ v; b v cd 



264 + -64^^ 



264 + -b^EW 



with 



.Gi 



— i — fBsUfn , B KR — i-j=f B fKfK ■ 



b's are the annihilation coefficients defined as 11 



h = %CiA\ , 

c 



C A A\ + C % A\ 



Cf 



C w A\ + C s Ai 



For the case of two pseudoscalars in the final states, A\ 2 are found to be 

'2- 



A\ ~ A\ ~ 2ixa s 9 



X A -A + 



7T 



+ r 2 x X 2 A 



(2) 



(3) 



(4) 



(5) 



where approximate SU(3) flavor symmetry has been used for r£ = r£ = r x ~ 2m K / {mb{m q + 
m s )) and Xa parameterizes the endpoint singularity as 



X A = ln 



m B 



(6) 



0.5GeV 

Notice that <$a is an arbitrary strong phase and normally pa < 1 is assumed, which reflects 
our ignorance on the annihilation amplitudes dominated by the soft gluon interaction. 
It is then straightforward to estimate the pure annihilation B decays in QCDF. Taking 



f Bs = 230 MeV 



m b (m b ) = 4.2 GeV 



m s (2 GeV) = 80 MeV 



(7) 
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FIG. 1. Contour plot of the branching ratio of B s — > tt + tt~ decay as a function of the annihilation 
parameters pa and 4>A- The solid red line represents the experimental central value and the light 
blue (grey) region corresponds to one sigma contour. 



and the Wolfenstein parameters 



A = 0.812 



A = 0.2254 



p = 0.144 



fj = 0.342 



(8) 



we show in Fig. [TJthe contour plot of B(B S — > 7r + 7r~) as a function of the parameters Pa 
and <pA- One may observe from Fig. [TJ that pa must be larger than 1.5 to account for 
the measured branching ratio of B s — > tt + tt~ within one sigma experimental error. Notice 
that we had chosen a relatively small s quark mass here which can enhance the penguin 
amplitudes, as adopted in the favored scenario S4 of [ll|. The large pa scenario has been 
discussed in [ll| and concluded to be unlikely as a fine-tuning of strong phase 4>a is required 
to satisfy the experimental bounds on ttK, ttK* and pK systems. But the experimental 
evidence of B s — > tt + tt~ forced us to have a closer look on the large annihilation scenario in 
the following. 

In the QCDF method, nonzero 4>a 111 the annihilation amplitudes could provide large 
strong phase which is also required to explain the measured direct CP violation of charmless 
B decays. Concerning the — > tt~K + channel which has been well measured 29] 



A CP (B d -> n-K + ) 



-9.8t\i)% 



B(B d -> n-K + ) = (19.4 ± 0.6) x 10" 



(9) 



one may determine, together with the B s — > 7r + 7r~ constraint, the annihilation parameters pa 
and 4>a as shown in Fig. |2j The form factor F BlT brings an important source of uncertainty 
to B d — > n~K + decay, which has been estimated in light-cone sum rules as 0.26 ± 0.03 
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FIG. 2. Contour plots of the branching ratio and direct CP violation of — > ir~K + channel, as 
well as the branching ratio of B s — > tt + tt~ , as functions of the annihilation parameters pa and <pA- 
The blue and green bands represent the branching ratio and direct CP violation of — > tt~K + 
decay, respectively, within two sigma experimental errors. The meaning of solid red line and the 
light blue region is the same as in Fig. [TJ The form factor F Bn is taken to be 0.26 (left plot) or 
0.22 (right plot). 



3fJ, |3l[- But as shown in the left plot of Fig. |2} there is no overlap between different bands 
if the central value of F Bn is taken. This means there is no solution of pa and <pA which 
can satisfy simultaneously the experimental constraints of — > n~K + and B s — > 7t + tt~ 
decays. Instead, a small (but still reasonable) form factor F Bn = 0.22 helps to reconcile the 
QCDF predictions with observations, as can be seen in the right plot of Fig. |2l One may 
also observe from Fig. [2] that, considering Bd — > ti~K + constraints, it is hard in QCDF to 
obtain B(B S — > 7r + 7r~) as large as the CDF central value 0.57 x 10~ 6 , though it is possible 
to reach the lower side of the 1 a error band. In the following we will take the parameters 



Pa 



1.75 , 



-53 c 



F Bn = 0.22 , 



0.28 , 



■^B S K 



0.26 , (10) 



where slightly small form factors F BK and F BsK have also been adopted, compared with 
the light-cone sum rules estimation 0.33 ± 0.04 jaOj and 0.30±g;g| fl, respectively. The 
superscript "PP" in the above equations means these annihilation parameters are adopted 
only for hadronic B decays into two light pseudoscalar mesons in the final states. For the 



Mode 


QCDF 


Experiment 


Mode 


QCDF 


Experiment 


B(B S -> 7T+7T-) 


0.40 


0.57 ±0.18 


B(B d 


0.20 


0.23 ±0.14 


7T-JC+) 


20.6 


19.4 ±0.6 




-11.4 


-9-81};? 


S(£+ -»• vr°^+) 


12.5 


12.9 ±0.6 


A C p(B+ -»■ 7T°K+) 


-2.3 


5.0 ±2.5 


-> 7r°ir°) 


9.4 


9.5 ±0.5 


6(5+ -»• vr+K°) 


23.2 


23.1 ± 1.0 




7.1 


5.0 ± 1.1 




37.4 


39 ± 17 


,8(5, -»• iT+K-) 


45.1 


26.5 ±4.4 


£(5 + -»• vr+vr ) 


5.8 


5.9 ±0.4 


iB(.B d -»• 7T+7T-) 


6.1 


5.2 ±0.2 


Acp(B d -)• 7T+7T") 


29.7 


38 ±6 


,8(5+ -»• K+if ) 


1.6 




£(5 d -»• JT if ) 


1.9 


96+ 21 

u - yD -0.19 



TABLE I. CP-averaged branching ratios (in unit of 10 6 ) and direct CP asymmetries (in units 
of 10 -2 ) of some B — y PP decay modes in the framework of QCDF, with the input parameters 
displayed in Eas. (l7l8ll0lll|) . 



parameters of wave functions which have less impacts on our results, we simply take 111. 133] 
A B = 200 MeV, 4 = 0.25, of = 0.06 , af = 0.25 . (11) 



We shall use the input parameters listed in Eqs.( 



as an illustration to check the 



QCDF predictions in comparison with data for some selected decay channels. 

Our results are listed in Table d for the CP-averaged branching ratios and direct CP 
asymmetries of some B — y PP decay channels. Notice that we have not discussed hadronic B 
decays with final states containing r\ or r]', as there are additional large uncertainties in QCDF 
relating to the relevant flavor-singlet components. For the decay mode B + — y tt°K + , naively 
one might expect Acp(tc~ K + ) ~ A. 



by 5.3cr. So we have followed 



16 



[n K + ) which is in disagreement with the experiments 



34| to adopt the scenario of large color-suppressed tree 
topology which may arise from spectator scattering or final state interactions. This Kit 
puzzle may also be explained by the so-called Pauli blocking effect proposed recently by 



Lipkin 



35| 



The QCDF results in Table |T] show good agreement with the experiments in general, 
except B s —y K + K~ and B d —y K°K° decays, whose branching ratios are estimated to 
be nearly twice larger than the experimental data. Notice that these two decay channels 
have something in common: both of them are penguin dominated with the annihilation 
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FIG. 3. Contour plots of the branching ratios of Bj — > tt~K + , B s — > K + K~ and B s — > 7r + 7r~ 
decays, as functions of the annihilation parameters pa and 4>a- The blue and green bands represent 
the branching ratios of B^ — > ir~K + and B s — >• K + K~ decays, respectively, within one sigma error. 
The meaning of solid red line and the light blue region is the same as in Fig. [TJ 

amplitudes determined essentially by a combination of 63 + 26 4 . For the other penguin- 
dominated modes listed in Table [H namely B — > irK channels, the annihilation amplitudes 
are largely determined by 63 term while pure annihilation decay B s — > 7r + 7r~ is dominated 
by 64 term. We show in Fig. [3] the contour plots of B(B d — > n~K + ), B(B S — > K + K~) and 
B(B S — > 7r + 7r~) in the plane of annihilation parameters pa — <Pa- One may observe that 
there is no overlap between different bands, which means there is no solution in the pa — <Pa 
plane, with other parameters fixed, to reproduce these three decay channels in agreement 
with the data simultaneously. For instance, one may take 

pY = 2 , 07 = "70° , (12) 

to get 

B(B S K+K~) = 29.6 x 1(T 6 , B(B S Tr+Tr") = 0.4 x 1(T 6 , (13) 
which are within one sigma error of the experimental data. But with the same parameters 
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we obtain 



B(B d tt-K + ) = 12.6 x 10~ 6 , A CP {B d -»■ 7r~A+) = -20.7% , (14) 

which are in strong contradiction with the experimental results (19.4 ± 0.6) x 10 -6 and 
—9.8t{'l%, respectively. 

We have tried to vary some sensitive parameters to find a solution. For instance, one 
may raise the form factor F BsK to have less tension with constraints of Fig. [3j but then 
B(B S —7- 7i + K~) will become too large. One may instead take a larger form factor F Bn but 
as we have shown in Fig. it can not work when AcpiBd tt~K + ) is included. We 

have also tried to vary the s quark mass from 80 MeV to 95 MeV but it does not help to 
reduce the discrepancy either. 

In QCDF calculations, we have followed the common practice to assume universal anni- 
hilation parameters pa and &a for all B ds — >■ PP decay modes, which is respected in the 
limit of SU(3) flavor symmetry. However, the SU(3) breaking effects could be as large as of 
0(20%). Observing that the difference of the annihilation parameters between Eq. tjTOl) and 
Eq. ffT2l is just about 20 percent, one possible way out of this problem is to introduce SU(3) 
breaking effects into annihilation parameters. That is to say, somewhat different pa and 4>a 
may be introduced for hadronic B decays into different final states. 

In the annihilation diagrams, the gluons may emit either from initial-state parton (de- 
noted as A\ 2 - i ) or from final-state parton (denoted as Ay. As a common practice, the 
annihilation parameters Xa has been assumed to be universal for both annihilation topolo- 
gies. But it is possible that Xa in is different from that in A\ 23 , as they originate from 
different topologies. Observing that B s — > 7r + 7r - depends only on A\ 23 , while B — > 7rK, 
K°K° and B s — > K + K~ decays contain both annihilation topologies but with different 
expressions, the above mentioned disagreement may also be solved by assuming X f A to be 
somewhat different from X\. As the SU(3) breaking effects of 0(20 — 30%) should be 
considered in general, it could be that A^'s are different for different decay channels and 
different annihilation topologies. 

In any case, as the QCDF predictions are very sensitive to the annihilation parameters, 
the predictive power in the framework of QCDF may be rather limited for many charmless 
decay channels. Notice also that, as there is no SU(3) flavor relation between B — > PP 
and B — > PV, VV decays, the large annihilation scenario in B — > PP decays does not 

8 



necessarily mean that pa should also be around 2 in B — > PV and VV decays. 

In summary, the first CDF evidence of B s — > tt + tt~ decays, which could be verified very 
soon by the LHCb collaboration, implies a large annihilation scenario with pA around 2 in 
the QCDF method. This is surprising as previous studies in QCDF preferred p& — 1- So we 
checked in details whether the large annihilation scenario is consistent with the experimental 
data for many well-measured B — > PP decay modes. Considering — > n~K + constraints, 
we observed that it is hard in QCDF to obtain B(B S — > 7r + 7r~) as large as the CDF central 
value 0.57 x 10~ 6 , though it is possible to reach the lower side of the 1 a error band. We 
found in addition that, taking slightly smaller form factors than the light cone sum rules 
estimation, the QCDF predictions are in good agreement with the data in general, except 
B s — > K + K~ and Bd — > K°K° decays whose branching ratios are predicted to be almost 
twice larger than the experimental measurements. One possible way to solve this problem 
is to take into account the SU(3) breaking effects in the annihilation parameters, which 
however means that the predictive power is rather limited in the QCDF method. 
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